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Abstract 

7-H ■ 

■ It is proved that small periodic singular perturbation of a cylindrical waveguide surface may 

On ' open a gap in the continuous spectrum of the Dirichlet problem for the Laplace operator. If the 

f**"* , perturbation period is long and the caverns in the cylinder are small, the gap certainly opens. 
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On 

O ■ 1 Spectra of cylindrical and periodic waveguides 

> : 

^ ■ 1.1 The cylindrical waveguide. 

Let ft = u> x K be a cylinder with the cross-section w C R 2 bounded by a simple closed contour dui 
assumed to be C°° - smooth for simplicity (cf. Remark [1] below). Interpreting ft, e. g., as an acoustic 
waveguide with the soft wall dfl, we consider the Dirichlet problem for the Hclmgoltz equation 

— A x v (x) = fj,v (x) , x G 0, v (x) =0,i£ dfl, (1.1) 

where A^, is the Laplacian, v the pressure, and /i a spectral parameter, proportional to square of the 
oscillation frequency. 

It is known (cf. [HH]) an d can be directly verified that, above a certain cut-off /z-f > 0, i. e. for 
/U > /Lif , the problem (jl.ip admits a solution in the form 

v (x) = exp {±i(z) V (y) (1.2) 

where i is the imaginary unit, z = x$ and y = (yi, 2/2) = X2) while 

M = (4 — C 2 (1.3) 
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and V arc an eigenvalue and the corresponding cigcnfunction of the model problem in the cross-section 

- A y V (x) = MV (y),y€u,V (y) = 0, y 6 dw. (1.4) 

Let Mi be the principal eigenvalue in the spectrum of the problem (|1.4j) : 

< Mi < M 2 < M 3 < ... < M fc ... -> +oo. (1.5) 

By the maximum principle (see, e.g., |30j). the eigenvalue Mi is simple and the cigcnfunction V\ can be 
fixed such that 

||Vi;L 2 M|| = 1, V 1 (y) > 0, y e u, d n V x (y) <0,y€ dw, (1.6) 
where d n stands for differentiation along the outward normal and L 2 (uj) for the Lebcsguc space. If 

M>/*t= M i> (1-7) 

then £ is a real number in ([1.311 . function (jl.2p does not grow or vanish as z — > ±oo and implies a wave 
which oscillates in the case ( ^ and stays constant in z for £ = 0. In other words, the wave propagation 
phenomenon occurs above the cut-off 

The problem (|1.1[) gives rise to the unbounded positive and self-adjoint operator Aq in L 2 (f2) with 
the differential expression — A x and the domain 

V (An) = H 2 (fi) n (fi; 90) . (1.8) 

We use the standard notation for the Sobolev space and the subspace H 1 (fl; 9Q) of functions in if 1 (0) 
satisfying the Dirichlct conditions in (jl.ip . 

The existence of the nontrivial wave (|1.2I) means that the point /i belongs to the continuous spectrum 
Cc (^.n) of the operator Aq. Indeed, multiplying v with the plateau function Xn (Figurc[T|) we see that 

\\X N v;L 2 {Sl)\\ 2 > 2{N - l)mes 2 w, || (A x + (j) (X N v) ; L 2 (O) || < const, NeN, (1.9) 

where N := {1, 2, ...} , and, therefore, i N~^XhV > is a singular Weyl sequence of An at the point u 

I > NGN 
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Figure 1: The plateau function 

whilst fx belongs to the essential spectrum a e (Aq) (see, e.g. [31 §9.1]). We emphasize that, by a general 
result in [18] (see also [T9l §3.1]), the kernel of the mapping 

H 2 (0) n i? 1 (0; 90) 9 w ^ - (A x + (jl) v e L 2 (Q) (1.10) 

stays finite-dimensional for any u € C and, hence, ct c (Aq) = a e (Aq) . 

The set C\ {/j, E C : Re /i >/-i|, Im/.i = 0} in the complex plane is the resolvent set of the operator 
Aq where the inhomogeneous Dirichlct problem 

- A x v (x) - fiv (x) = f (x) , v (x) = 0, x e dn, (1.11) 
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has a unique solution v G H 2 (fl) for any right-hand side / G L 2 (f2) and the attendant estimate 



\v;H 2 (Q)\\ < Cp l|/;L 2 (f7)| 



(1.12) 



ensures that the mapping (jl.lOp is an isomorphism. In contrast, on the continuous spectrum, mapping 
(ll.lOP looses even the Fredholm property (cf. [T9l Thm. 3.11]) and the inhomogeneous problem (| 1 . 1 1 [) 
requires a specific formulation involving radiation conditions at infinity. In the sequel we do not need 
such the formulation and refer, e.g. to [HE], [HI §5.3] for details. 



1.2 The periodic waveguide, a quasi-cylinder. 

Let II be a domain in R 3 with a periodic cross-section (Figure [2]) . More precisely, II is the interior of 
the union 

n = U^> a-13) 

where Z= {0, ±1, ±2, ...} , Wj — {x = (?/, z) : {y,z — j) G tu} , and the reference periodicity cell w lies 




Figure 2: Periodic waveguide 



inside the circular cylinder {x : \y\ < R, z G (—5, 5)} of radius R > and the unit height. We, of 
course, assume that II is a connected set, i.e., a domain. 

Due to possible boundary irregularities the Dirichlet problem for the Helmgoltz equation in the 
quasi-cylinder II needs the variational formulation as the integral identity [T7] 



(V x u, V x v) n = A (u, «) n , v g H 1 (H; an) . 



(1.14) 



where X7 X is the gradient operator, ( , ) n the natural inner product in L 2 (II), and A a spectral parameter. 
If the boundary <9II is smooth, the integral identity (|1.14[) with arbitrary test function v G H 1 (II; 911) 
is equivalent to the differential problem of type (|l.lj) in II. The spectral problem reads: To find fj, G C 
and a nontrivial function u G H 1 (0; d£l) verifying (|1.14| . 

Remark 1 We could take any open connected and bounded cross-section uj of the cylinder U, and for- 
mulate a spectral variational problem of type \ 1. 1$ in Q. However, the smoothness of 8lo will be used 
in § 3 for an asymptotic analysis. 



1.3 The band-gap structure of the essential spectrum in a quasi-cylinder. 

The left-hand side of (|1.14[) is a positive continuous form in H 1 (II; dH). According to [31 §10.2], this 
form is associated with a positive self-adjoint unbounded operator An in L 2 (II). If the surface 911 is 
smooth, An gets the same properties as Aq with the only exception, namely, its essential spectrurrf^ has 
the band-gap structure 

^n) = U peN T P (1-15) 

where T p are closed segments 

T P =[A;,A+]. (1.16) 

1 The authors do not know if it is possible that a segment T p collapses into the single point Ap = Ap which thus 
becomes an eigenvalue of the operator Api- In the case Ap > Ap for any p = 1,2,... the essential spectrum cr e (Ari) 
coincides with the continuous spectrum cr c (An) as in the straight cylinder. 
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Formulas (jl.ISp and ()1.16j) remain valid without the smoothness assumption (see [10j[T5] and others). 

To indicate the segments (|1.16[) , the model spectral problem on the periodicity cell vj must be 
considered 

Q v (U, V; vj) := {{V x + ir,e 3 ) U, (V, + i V e 3 ) V) w = A (U, V)„ , V G H x per {vj; 7) , (1.17) 

where ej is the unit vector of the 2 j -axis and H per {vj; 7) is the subspace of 1-periodic in z functions 
V G H 1 (w) vanishing on the lateral side 7 = { x G dw : z G (— \ , |) } of the cell. If, for certain 77 G [0, 2ir) 
and A > 0, the model problem (|1.17[) has a nontrivial solution U G H per (0757), then the Flochet wave 

u (y, z) = exp (irjz) U (y, z) (1.18) 

satisfies formally the original problem in the quasi- cylinder II that is the integral identity (|1 . 14[) with 
any test function v G (n) (infinitely diffcrentiable functions with compact supports). One readily 
constructs from the Flochet wave the singular Weyl sequence for the operator An at the point A = A 
with the help of the plateau function drawn in Figure [1] (cf. formulas (|1.9p ). 

For any real 77, the sesquilinear form on the left of Q1.17P is Hermitian. closed and positive. Thus, 
problem (|1.17p can be associated with the unbounded self-adjoint positive operator .An (r)) in L 2 {vj) 
(see again [3j §10.1]). The domain V {An (??)) is included into the Sobolev space H 1 {vj) and, therefore, 
is compactly embedded into L 2 {vj). By [31 Thm. 10.1.5], the spectrum of _4n (77) is discrete and forms 
the infinitely large sequence 

< Ai (77) < A 2 {rf) < ... < A p (77) < ... -> +00 (1.19) 

of eigenvalues which are listed according to multiplicity. The functions 77 — > A p (77) are continuous (see 
|27[ Ch. 9]) and, by an evident argument, 2ir - periodic. This means that the endpoints of segments in 
()1.16p are calculated as follows: 

A± = ±max{±A p (r7) : 77 G [0,2tt)}. (1.20) 



1.4 The Fourier and Gel'fand transforms. 

Let us comment on the above-mentioned inference. A correspondence between the problems (|1.1[) in 
the cylinder 51 and f)l .4() in the cross-section u> is pointed by the Fourier transform (see [18] and e.g. 
[19\ I20j). For the quasi-cylinder II, one ought to apply the discrete Fourier transform, namely, the 
Gel'fand transform 

u (y, z) — > u (y, z; 77) = -j=^ cxp (-777 {z + j))u {y, z + j) (1.21) 



2tt . c 

(see [2H] and e.g. [HI US]). Note that (y,z) G II on the left of (TOT]) but (y,z) G vj on the right. The 
Gel'fand transform establishes the isomorphisms 

L 2 (n) » L 2 (0, 2tt; L 2 {vj)) , H l (n) » L 2 (0, 2n; H l per {vj)) 

where L 2 (0, 2ir; 03) stands for the Lebesgue space of abstract functions, 



( > 
J ||C/(77;23)|| 2 d77 



and 03 is a Banach space. The corresponding Parceval theorem provides the identity 

2ir 



u{x)v{x)dx = / u{x;rj)v{x;r])dxd'r], u,v G L (II) 
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which, together with the formulas 

v(y,0;r]) = v(y, l;r/), d z v(y,z;r)) = d z v(y,z;r])-ir)v(y,z;ri), v £ C™° (IT) , 

indicate the immediate correspondence between the problem (|1.14p in the quasi-cylinder II and the 
family (77 £ [0, 2n)) of problems ()1.17|) in the periodicity cell w. 

A result in [TB] (see also [101 [T51 [TJ5] and others) demonstrates that the operator of problem (|1.14[> 
with the fixed A S C regarding as the mapping 

h 1 (n ; an) — ► h 1 (n ; an)* (1.22) 

is Frcdholm if and only if, for any r\ £ [0, 27r), the problem 

((V x + i??e3 ) (7, (V, + i V e 3 ) V)„ - A (U, V)„ = T (V) , V £ (w; 7) , 

with the fixed A G C is uniquely solvable with any linear functional T £H^ er (zu\ 7)* on the space 
Hp er (m;-f). The fact mentioned above ensure the segmental structure (| 1 . 1 5|) of the essential spectrum 
a c (A n ) and formulas (fi~16|) . (TOO)) , (fTTT9)l for the segments. 

1.5 Gaps in the essential spectrum. 

The band structure (| 1 . 1 5[) of the essential spectrum a e (An) in the periodic waveguide n allows for gaps, 
i.e., intervals on the real positive semi-axis R+ which lie outside a e (An) but have both the endpoints in 
<r e (An) ■ As was commented, such a gap cannot appear in the essential spectrum a e (Aa) = a c (An) of 
the cylindrical waveguide VL. However, the segments (|1.16|) for the periodic waveguide can intersect each 
other and, as a result, cover the whole ray [A f , +00) . In other words, even a quasi-cylinder can have the 
essential spectrum with the only cut-off A t and no gap. 

The main aim of the paper is to show that a small periodic surface perturbation of the cylinder 
O = w x K opens a gap in the essential spectrum of the corresponding perturbed quasi-cylinder H h . 

In the literature results on opening gaps are mainly related to periodic media in the whole space 
IR m of a piece-wise constant structure described by either scalar differential equation, or the Maxwell 
system. We refer to papers [H [5j [6j [71 [9] and reviews [TUl [TTj . Usually the existence of a gap in the 
essential spectrum is established by assuming contrast properties of the media and selecting or matching 
the coefficient constants. Results of a different kind are obtained in p21 H31 E] and the present paper, 
namely coefficients of differential operators are constant and invariable but gaps are opened by varying 
the shape of the periodicity cell forming a quasicylinder. In [121 1 1 3j two-dimensional periodic waveguides 
of thin width are investigated while spatial waveguides with either regular perturbation of a cylindrical 
boundary, or a periodic nucleation are studied in [14j . 

2 Opening a gap in the continuous spectrum in the perturbed 
periodic waveguide 

2.1 Any cylinder is a periodic set. 

The straight cylinder fl = uj x R can be regarded as the quasi-cylinder n° with the periodicity cell 
w° = uj x (0, 1) . The wave v (x) = exp (i(,z) V (y) (see (|1.2|) ) turns into the Floquet wave (|1.18|) with 
the attributes 

rj = C - 2jrq (() , U (y, z) = exp (2mq (Q z) V (y) , (2.1) 

where q (() = m&x{q £ Z :2irq < £}. We point out that the factor exp (2niq (() z) is 27r-periodic in z. 
Thus, each of the curves 

,i = M p -C 2 , (2.2) 
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forming the continuous spectrum a c (fi) (cf. (| 1 .3|) and ()1.7j) ). gives rise to infinite number of pieces 

A = M p + (?/ - 271-g) 2 , gGZ,ije[0, 2tt) (2.3) 

generating segments in ()1.15|) which cover the whole ray [Mi, +00) as it is shown on Figure [3] for the 
lowest curve (|2.2p with p = 1. In particular, under the assumption 

Mi + 7T 2 < M 2 (2.4) 

(see Remark [2] below) the spectral problem (JTTTTJ) in the cylindrical cell tu°=cjx(— \,\) gets the first 



Figure 3: Constructing Floquet waves in the straight cylinder. 



cigenpairs 



a? fo) 



Mi + r/ 2 



»? e [0, tt) , 



Mi + (2tt - , £ (tt, 2tt) . 



C/ 1 °(y,z; 7? ) = y 1 (y) 



1, 



?7 e [o,7r), 



exp (—2-Kiz) , 77 e (7r, 27r) , 



while, at 77 = 7r, the eigenvalue 



A? (tt) = A" (tt) = M a + tt 2 
becomes of multiplicity 2 and has the cigenfunctions 

U° (y, z) = Vi (y) , C/° (y, z) = Vi (y) exp {-2mz) . 



(2.5) 
(2.6) 

(2.7) 
(2.8) 



It is known (see, e.g., [2T1 §7.6] and [2T1 Ch.9.10]) that a small perturbation of the cell w° prompts 
perturbations of eigenvalues in p,19[) . Two situations drawn in Figure 3] may occur for the first couple 
of eigenvalues and we shall show that a periodic singular perturbation of the cylinder f2 (see Figure [5]) 
provides opening a gap in the continuous spectrum (as indicated by over-shadowing in Figure [He). 





Figure 4: The perturbations (b and c) of the eigenvalue curves (a). 
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Figure 5: Singular perturbation of the periodicity cell. 



2.2 The singular perturbation of the cylindrical surface. 

To describe the boundary perturbation of the cylinder II = fl, we introduce in a neighborhood T of 
the contour T = dcu the natural curvilinear coordinate system (n, s) (Figure [6]) where n is the oriented 
distance to T, n < inside wflT and s is the arc length on T evaluated from a point O' € T counter- 
clockwise so that the point O = (C, 0) € dH° has the coordinates n = 0, s = 0, z = 0. 



Figure 6: The curvilinear coordinates. 



Given a small parameter h, we introduce the sets 

6 h = {x e T x R: £ :=hr x (n, s, z) e 0} , w' 1 = zu°\¥ (2.9) 

where 6* is a bounded nonempty domain in the half-space R?_ = {£ = (£1,^2, £3) : £1 < 0} . According to 
formula (|1.13p , the reference cell w h in (|2.9p generates the quasi-cylindcr n' 1 with a singular perturbation 
by the 1 - periodic family of the caves or superficial voids 9- = {ielxl: ft, -1 (n, s, z — j) 6 0} (Figure 
0a,b). 

Remark 2 IFe /iGwe assumed that the perturbation period T is equal to 1. J/ T / 1, the rescaling 
x 1 — > T~ 1 x turns the cylinder fi into Qt = wj- x R w/iiZe i/ie model problem <\l-4ty in the new cross- 
section lot = {y ■ Ty € to} gets the eigenvalues T 2 Mk where Mk are taken from (j 1 . 5} . Since M\ < M2, 
the assumption <\2.4l is satisfied in the case 

T > 7T (M 2 - Mi)"' . (2.10) 
2.3 The boundary layer phenomenon. 

To examine the behavior of eigenfunctions in the periodicity cell w h near the boundary perturbation, 
we need to construct the boundary layer (see, e.g., [35], [HJ Ch. 2.9]). To this end, we use the stretched 
coordinates £ in (|2.9[) . Since the Laplacian A x in the curvilinear coordinates reads 

A x = (1 + nxOO)- 1 (J- (1 + n* (-)) £ + + „x (s))" 1 1) + (2.11) 
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where x(s) is the curvature of T at the point s, we formally have 



A x ~ h~ 2 A i + nip') — - 2h(0')£x-pz 



(2.12) 



Hence, in view of formulae (|2.9p , the coordinate dilation x 1 — ► £ leads to the following limit problem 

- a s w (0 = 0, f g e, tu (0 = .9 (0 , £ e ae, (2.13) 

in the imperfect half-space (Figure 0a, b) 

e = K'i\£. (2.14) 
It is known that, for a sufficiently smooth datum g with a compact support, problem t|2.13p has a unique 




Figure 7: The rescaled perturbation of the boundary. 



solution with a finite Dirichlct integral. In the sequel we need such the decaying solution W (£) with the 
special right-hand side g (£) = — £1 which vanishes on 80\d9 and obeys the asymptotic form 



where R > is fixed such that Ixl < R for x £ 9. 



(2.15) 



Note that — ^2-7r |£| J £j implies the Poisson kernel and Pg > by virtue of the maximum principle. 

Remark 3 The exterior Dirichlet problem for the symmetrized set 9" = |£ : (— |£i| ,£2,^3) £ ^} ( c f- 
Figures^ and\$j) has an intrinsic integral characteristics, the polarization matrix (see \29l Appendix G]), 
which is extracted from asymptotics at the infinity of the harmonics Wj under the Dirichlet conditions 
Wj (£) = — £j, £ € 89". The odd extension ofW from onto R 3 \#" coincides with W3 and, therefore, 
Pe is proportional to an entry in the polarization tensor of 9" . We call Pg the polarization coefficient 
of the cavity or void 9 in the half-space. 




Figure 8: The symmctrization of the boundary perturbation. 
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2.4 The main result on asymptotics. 

To identify the gap, we need two assertions on eigenvalues and cigenfunctions of the auxiliary problem 

((V x +irie 3 )U h ,(V x +i V e 3 )V)^ h =A h {U h ,V)^ h , V £ H^ r {zu h ; j h ) , (2.16) 

in the perturbed periodicity cell vo h in (|2.9[) with the lateral side ^ h = {x £ dw h : \z\ < 1/2} . We 
enumerate the eigenvalues in the same way as in (| 1 . 19[) : 

< A? (rj) < A h 2 (r,) < ... < A p (r?) < ... - +co. (2.17) 

However, under the assumption (|2.4|) the first couple of eigenvalues in (|2.17|) is denoted by A^_ (n) while, 
according to (|2.5[) . we have A^_ = Mi + (rj — it ± 7r) 2 in the limit (/i = 0) problem (|1.17|) in zu = m° and 
the corresponding cigenfunctions are given by (|2.6|) . 



Theorem 4 There exist positive numbers ho,(3o,CQ such that, for any h £ (0, ho] and \[3\ < [3oh~*, the 
first couple of eigenvalues in <\2. 11\ of the problem l\2.16} on the periodicity cell w h , determined in 
takes the asymptotic form 



A'l (tt + (3h 3 ) = Mi + 7T 2 + h 3 [V±^P 2 + iTT 2 l3 2 j +Ai(ir + f3h 3 ) (2.18) 
where the remainder admits the estimate 

Aj_ (ir + ph 3 ) < C A hi (2.19) 

and the positive quantity 

V =P e IdM (O'f (2.20) 

is calculated according to \2. 15\ and 

The asymptotic formula (|2.18|) . (|2.20p will be derived in §3 and the remainder estimate (|2.19p in §4. 
To detect a gap in the continuous spectrum of the problem (|1.14[) in the quasi-cylindcr II , we also prove 
the following intelligible inequalities. 

Lemma 5 Entries of the eigenvalue sequences <\2.17^ and (| of the auxiliary problems in the cells 

w h and w° , respectively, are in the relationship 

A° p (v)<A h p ( v )<A°(n)+C p h 3 , (2.21) 

where C p is independent of rj £ [0, 27r) and h £ (0, ho] ■ 

Proof. Let A^ be a unbounded operator in L 2 (za h ) generated by the closed positive Hermitian 
form Q v (•,-,tu' 1 ) on the left of (|1.17[) (cf. [3l §10.2]). We employ the max-min principle (see [3l Thm. 
10.2.2]) 

A^ (rj) = max inf ^" V 1, p £ N. (2.22) 

P £p ue£ P \{o}\\U;L 2 (m h )\\ 2 

Here £ p is any subspace in Hl er {■co h ; r y h ) of co-dimension p — 1, in particular, £\ = H per (m h ;^ h ) . 

Let the eigenfunctions U p (-,rj) corresponding to A p (rj) satisfy the normalization and orthogonality 
conditions 

([/;,C/J) Bl =i M)M eN, (2.23) 

where S p , q stands for Kronecker's symbol. The subspace 3 P C H per (vu h ; 7' 1 ) is spanned over the functions 
X h U?, ...,X h U° while X h £ C™ r is such that 

X h (x) = for \x-0\< C x h, X h (x) = 1 for \x - 0\ > 2c x h, (2.24) 
X h (x) = for x£ 9 h , 0<X h < 1, \V x X h (x)\ < ch' 1 . 
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In other words. Xh is equal to 1 everywhere in w h , except in the vicinity of O, and Xh vanishes in the 
cavern. We have 

{X h lP, X h U*)„ h = {UP, U*)^ + ((1 - Xl) U", U<>) w0 > S p , q - c pq h 2 h 3 , (2.25) 
Q n {UP, U"; w h ) < Q v {IP, U"; w°) + c pq (\\U*; H 1 {E h )\\ + ftr 1 \\lP- L 2 {E h )\\) • 

■{p^H 1 (S,)||+/.- 1 ||C/P;i 2 (S h )||) 
<A° p 6 p , q + c pq {h 3 + h~ 2 h 2 h 3 ) . 

Here come the factors ft -1 and h 2 from the differentiation of Xh and the formula 

\U P {x)\ 2 < Cp \x- Of < C p h 2 , x G E h , 

while h 3 is order of the volume of the set Eh = supp(l — Xh) D supp|V x ^h| • 

The intersection of the subspaces £ p and J p contains the nontrivial linear combination 

p p 

U p = X h J2a P U p , J2 K\ 2 = L ( 2 - 26 ) 

j=i j'=i 

Hence, according to (|2.22p and (|2.23[) . we derive that 



p 



Ap (rj) < max — < 



]>>° (r?) + C p h 3 



£ P \\U p ;L*{m h )\\ 2 ~ 1-C p h 3 

and the right inequality in (|2.21[) is proved. 

The left inequality can be easily derived by applying the max-min principle to the operator A® and 
extending eigenfunctions Up by zero from w h onto zu°. 

2.5 Detecting the gap. 

By Lemma [5] and formula (|2.3p . we conclude that 

Ai l (r?) < Mi + min {r/ 2 , (2tt - r/) 2 } + Cift 3 , (2.27) 
A£ (??) > min |m x + max {t/ 2 , (2tt - t/) 2 } , M 2 + min jr? 2 , (2?r - ??) 2 }} ■ 

Hence, in view of the assumption (|2 .4(1 we can choose ijo > (27r) 1 max{Ci, 3T 5 } and fto > such that, 
for ft, G (0, fto] and 77 G [0, 7r — 770ft 3 ) U {it + f]oh 3 , 2ir) , the interval 

(Mi + 7T 2 , Ml + 7T 2 + 2-pft 3 ) 

is free of eigenvalues (|2.17p . In the case 1 77 — 7r| < 7/0 ft 3 we apply Theorem [4] to observe that also the 
interval 

' Mi + 7T 2 + C A hi , Mi + 7T 2 + 2-pft 3 - C A ft5 ) (2.28) 



does not contain the eigenvalues. We emphasize that the endpoints in (|2.28j) arc established by the 
following inequalities taken from Theorem 3) 



A1 1 (7/) < A/i + 7r 2 ft 3 \ V - Vyjl + 4ir 2 V- 2 v%j + C A h? < Mj + tt 2 + C A h*, 

A2 {v) > M i + ^h 3 \ V + V^Jl + A^V- 2 ^ - C A hi > Ah + n 2 + 2Vh 3 - C A hh. 



These two facts provide the main result in the paper. 
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Theorem 6 Under the assumption §2.4h there exist positive numbers ho and cq such that, for h € 
(0, ho], the essential spectrum (| 1.15\ of the problem Q 1. H\ in the periodic waveguide H with the period- 
icity cell zu h in <\2.9§ has a gap of length I (h), 

\l(h) -2Vh 3 \ < c hi, (2.29) 

situated just after the first segment in 15§ . Here V is the positive quantity $2.20\ . 

We finally mention that under the assumption Mi + n 2 > M2 opposite to (|2.4[) , the first and second 
segment and T§ intersect (cf. Figure [9] where two dotted curves correspond to (|2 . 3[) with p — 2 and 
q = 0, 1) and therefore, the gap discovered in Theorem [6] does not occur. This conclusion readily follows 

from the rough estimate (|2.2ip for the pcrturbated eigenvalue M2 +min |r/ 2 , (2n — r/) 2 | in the auxiliary 

problem in w . 




Figure 9: Overlapping of the bands. 

If M 2 = M 1 +it 2 , then the gap is still open because Aj> (0) > M2 and Theorem [4] is still valid. At 
the same time, the gap length is O (ft 3 ) only in the case when M 2 is simple and d n V2 (C) 7^ 0, but 
I (h) = o (/i 3 ) provide the normal derivative of an eigenfunction corresponding to A/2 vanishes at the 
point O' . This conclusion can be confirmed by an asymptotic analysis of eigenvalues, similar to §3 and 
§4. However, calculations become much more combersome and we omit them here while refereing to \2\\ 
Ch. 9,10] for general asymptotic procedures. 

3 The asymptotic analysis 
3.1 The asymptotic ansatze 

Let us examine the eigenvalues A± (77) of the spectral problem (|1.17p in the perturbed periodicity cell 
w h which are close to the double eigenvalue (|2.7|) of the problem in zu°. We fix the dual variable of the 
Gel'fand transform 

77 = 7T + /% 3 , (3.1) 

where (3 £ R is the deviation parameter. By varying (3, we watch over the eigenvalues A^_ (77) in the 
vicinity of the collision point in Figure [3Ja. Note that the factor h 3 is adjusted with the second term in 
the eigenvalue asymptotic ansatze [23] (see also [U Ch. 9] and [55]) 

Al(r))=A + h 3 A' ± (f3)+A>t(r,). (3.2) 

Here A = M\ + ir 2 , A ± (/?) is a correction term to be found out and A^. (77) a small remainder to be 
estimated in The asymptotic ansatze for the corresponding cigenfunctions looks as follows: 

Ui (x; n) = Ul (x; (3) + h x (x) («4 (£; (3) + hw 2 ± /?)) + h 3 U' ± (x; p) + U h ± (x; 77) . (3.3) 
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The main term 

U° (x; (3) = a\ ((3) U+ (x) + (/?) U- (ar) , a = +, - (3.4) 

is a linear combination of functions (|2.8[) with the coefficient column a 0- = (a? , al) T while |a CT | = 1 and 
T stands for transposition. The boundary layer terms w% are intended to compensate for a discrepancy 
produced in the Dirichlet condition on the surface d9 h n dw h by the term U±. Since w a (t;;{3) is 
defined only in the set U n , the cut-off function \ is introduced in (|3.3p such that x = outside 
a neighborhood of O and x — 1 m the vicinity of the point O. The correction term JJ' a is used to 
compensate for discrepancy of U® and hxw\ in the equation with the differential operator 

A y + (d z + irjf + K h a (?]) ~A y + (d, + ivf + A + h 3 (2i0 (d z + iit) + A a (/?)) + ... (3.5) 

=:L°{V x ) + h 3 L / {V x ;P) + ... 

which is decomposed in accordance with (|3.1|) and (|3.2|) . Notice that the second boundary layer term 
hxw± is linear in jt {€)') (cf. Q2.12p ). however it does not influences A ± (/?) in (|3.2p and becomes important 
only in S]4] for justification estimates (see Section [3j4) . This observation displays the effect of opening 
the gap to be independent of the curvature h of the contour duj. 

The function [7 ± in (|3.4|) gets a singularity at the point O and, hence, we ought to introduce another 
cut-off function Xh into (|3.3|) (see (|4.8jl ). However, since the asymptotic analysis in this section is formal, 
we avoid to multiply U± with Xh here. To accept this mathematical licence, one can assume that the 
coordinate origin £ = lies inside 9, i.e. O ^w h for any h G (0, ho] . 



3.2 Calculating the asymptotic terms 



Since the eigenfunction V\ of the problem (| 1 .4|) is smooth near the boundary du>, the Taylor formula and 
the Dirichlet condition yield 



U° (x; (3) = {a% {(3) + a a _ (/?)) nd n V x {&) + J2 a ± W u± (»> ». *) + O (|n| 3 + |s| 3 J = (3.6) 

= ft « 09) + (/?)) &d n V (O') + h 2 Y,a a ± ((3) U± (£) + O (h 3 ) , 

± 

x e <96> h n cW\ 

(0 = f-^i (C) + fifaSa^i (C) - 2-Kt8 ± ^ 3 d n V l (C) . (3.7) 

Here we used the definitions of U® and £ in (|3.4p and (|2.8p . Recalling the special solution W of the limit 
problem (|2. 13[) with g (£) = — £i, we set 

u£ (0 = (a^ 09) + <£ (/?)) <?„Vi (O') W (£) (3.8) 

in order to compensate for the main discrepancy O (h) in (|3.6p . By the asymptotic expansion (|2.15p we 
obtain 

hut (0 = -M. (/?) ^3 + O QjS J = (0 ^ + O (£) (3.9) 

where r = V n 2 + ,s 2 + z 2 = h |£| and 

A„ {(3) = {a° + 09) + al (/?)) d„Vi (C) P e . (3.10) 

After applying the differential operator (|3.5p to the right-hand side of ()3.3[) we collect coefficients on h 3 
and derive the differential equation 

L {V x )U' rT {x ] (3)=F' (T {x;(3)~ (3.11) 
:= —L' (V*; /?) /?) + L ° (V«) (x (*) 4, 03) ^3) ,i6b , 
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which is to be supplied with the following Dirichlct condition on the lateral side 7 of the cylindrical cell 
w° : 

EC (*;/?) = 0, xe 7 °. (3.12) 

The first term L'U® on the right of ([3. lip is smooth in zu° but the second one gets a singularity at the 
point O £ 7 . By means of (|2.1ip (see also (|2.12p and (|3.5[) ) we conclude the representation 

where £° is a first-order differential operator. Hence, 

L° (V«) (x (x) r~ 3 n) = O (r" 3 ) , r -> +0. (3.14) 

The strong singularity (|3.14|) of the right-hand side does not allow for a solution U' a of problem p. lip , 
(f3~T2|) in the Sobolev space H 1 (n7°; 7 °) . 

3.3 The regular correction term for the eigenfunctions. 

Let us move into the scale of Kondratiev spaces V\ (see [18] and, e.g., [El [20]) equipped with the 

weighted norm 

\\U;Vl{™°)\\= (j2\\p T - l - k ^UlL^^)\\ 2 Y (3-15) 
\fc=o / 

where p(x) =dist(a;,0), V£?7 is the family of all order k derivatives of U while I £ {0, 1, ..} and reR 
are the smoothness and weight indices, respectively. By the one-dimensional Hardy inequality with the 
particular exponent a = 1, 

+00 +00 „ 

4 



p «-i\u(p)\ 2 dp<^ / p a+1 





we obtain that 



du ^ ^ 
dp 



dp, a > 0, u £ C\ [0, +00) , (3.16) 

\\p- x U- L 2 \\*<c\\U; H 1 (w°) f . (3.17) 
Thus, the space 

Vo, P er ( ro °> 7°) = {U £ (w°) : U = on 7 , U is 1 - periodic in z) 
coincides with the space H per (ci7 ,7 ) algebraically and topologically. This means that the mapping 

^(^7°)^oU(^7°)% (3.18) 
associated with the problem (|3.1ip , (|3.12p , inherits all properties of the mapping 

^ 1 er (^ ,7°)->i? p 1 e ,(^ ,7 )*. (3.19) 

Moreover, theorems in [TB] on lifting smoothness and shifting the weight indices (see also [THl Theorems 
4.1.2 and 4.2.1]) convey these properties to the mapping 

KU (-°, 7°) n V^ per «) -> V^ per «) (3.20) 

in the case 

'eN, re (-§,§). (3.21) 
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Remark 7 The bound ±| for the weight index r in <\3.21$ may be computed as follows: the "linear" 
function x i — > x( x ) n belongs to V^li. 1 under the restriction r > — | w/iiZe i/ie Poisson kernel 

x(x)nr~ 3 lives outside V^i 1 (^°) Ml £fte case r < |. explanation of such a mnemonic rule, main- 
tained by the general theory, can be found in the introductory chapters of books U!A \20§ . 

Taking into account the singularity of at the point O, we see that 

f^CM for any re 

Recall that A = Mi + tt 2 is a double eigenvalue of problem (| 1 . 1 7[) in the cylindrical cell zu° (see (|2 . T[) 
and (|2.8p ). Thus, the co-kernel of the mapping (|3.20p is spanned over the functions U± and the problem 
(|3~TTj) . ([3TT2]) admits a solution in (ru°) with re (5, |) if and only if 



t/± (x)i^ (x;P)dx = 0. (3.22) 

j 

Note that t/± (x) = O (|n|) in T x (— |, |) and, in view of (|3 . 14[) . the integral in (|3.22[) is convergent. 
Let the compatibility conditions (|3.22p be satisfied. The orthogonality conditions 



U± {x)U' a [x;(3)dx = (3.23) 



make the solution unique. 



Remark 8 . The general results \18[ \23\j , (see also U91 Ch.2.3]) furnish an asymptotic form of the 
solution U' a . By (fOIj) . $2~T2} and (EM, (Ell, we have 

F' a {x\(3) =r- 5 Q F (n,s,z) + 0(r- 2 ) , r -»• 0+, (3.24) 

where Q ... stands for a homogeneous polynomial of degree 2. A routine and traditional calculation brings 
the expansion 

U' a (x; f3) = r- 3 Qu (n, s,z) + (r°) , r ^ 0+ , (3.25) 

which as well as §3. 2$ can be differentiated under the convention V x O (r*) = O (r* _1 ). We need not 
explicit formulas for Qp and Qjj , however the estimate 

V k x U' a (x; (3) <c k (l + \[3\) r- 1 -", k = 0, 1..., (3.26) 

inherited from (| 3. 25\ . will be useful in Constants c k in (| 3.26$ do not depend on the parameter (3 
while Qp and Qu are linear in (3 (see formula for L' in (|ff.<5[) ). 

All the above conclusions, of course, are known explicitly for the Poisson kernel. 
3.4 The correction term for the eigenvalues. 

Let us compute the left-hand side of (|3.22p . Applying formulae (|3.5|) . (|1.6| and (|2.8| . (|3.4|) . we readily 
get 

h = f W{x)L' (V*; (3) U° a (x; (3) dx = (3.27) 

V(y)\ 2 dyJ (exp ((-tt ± tt) iz) (a+ ((3) (-2tt/3 + A' a (/?))) + 

+ (a- (/3) (2ir(3 + K' a ((3)) exp (-27riz))) dz = 
at (/?) (T2^/3 + A; (/?)). 
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To calculate the second integral, we employ the method [22]. Using the Green formula in the domain 
txj°\Bs where Bg = {x G T X (—5, 3) : r < <5} and <5 > is small, we have 

h= [ U^L°(V x )( X (x)A a -^)dx = A a lim [ U^L° (\7 x ) *j^dx = (3.28) 



<5-»0 



d n n d 



= -A CT lim/ [/ ±(a; ) _ t/ ±(a; ) dSx 



d d z 

Here — — = h «7r— and AT is the interior normal on the surface dB$nvj® . Since the gradient operator 

QN dN r 
in the curvilinear coordinates takes the form 

_,(l + nx( S )) -,-), 
we obtain 

— - 

N(x)= (r 2 +s 2 ((l + nx(s)r 2 -l)) 2 (n, (1 + (s))" 1 s, z) , 

A . + ((i + „(.))-> - 1))" 1 („A, (i + „ W)-4,.|) ■ 

Thus, computing the limit in (|3.28|) . we can make the changes 

TT , \ ^ <9 ^ S ^ Z ^ 

U± (x) 1— > n&Vi C ,^^7T = -ir + -^- + "TT- 

(cf. (|2.8j) . (|3.6j) ). Taking the relation mes 2 (<9£\5 n z&°) = 2irS 2 + O (<5 3 ) into account, we then arrive at 
the formula 

h = -A„(p)d n V 1 (0')lim f ( n l dSx = (3.29) 

= A a {13) d n V x {&) = (a+ (/?) + a- {(3)) Pg \d n V x {0')f . 

Here we have used notation (|3.10p and further we set V —Pg \d n V\ {0')\ 2 as in (|2.20p . 

By (|3.27[) and (|3.29|) . the compatibility conditions (|3.22[) reduce to the system of two algebraic 
equations 

±27r/3a± 09) + V {at (/?) + a~ (/?)) = K (P) 4 {(3) . 
Eigenvalues of the corresponding matrix 



2tt/3 + V V 

V -2TTf3 + V 



(3.30) 



look as follows 



A' ± {(3) = V±VV 2 + An 2 p 2 . (3.31) 



3.5 The second term in the boundary layer. 

Even in the case x{0') = 0, e.g., the contour dui is flat near the point O' and, by formulas (|1.1|) and 

(EH}, 

dlV^O') = -x(O')d n V 1 (O') = 0, 
the second term (|3.7[) of the discrepancy (|3.6[) does not vanish. The Sobolev norm of the functions 
h q x w ± i s O (h p+ ^ and, therefore, our aim to derive estimates with the bound chi forces us to deal 
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with h 2 xw± in §4, although this term, owing to the proper decay as |£| — > oo, does not influence the 
correction term h 3 A' ± ((3) in (pT2"|) . 

The Taylor formula (|3.6[) gives immediately the boundary condition 



w 2 ± (£; P) = -at 09) U+ (0 - a ± (/3) IT (0 , £ G 96. (3.32) 
To derive the differential equation 

x (OO — - 2x (O') Ci ^| + 2m— J wi (f ; /?) , £ G 6, (3.33) 

requires much more elaborated analysis based on the procedure [2lT §2.2. Ch.4] of discrepancies rear- 
rangement. First, the differential operator C\ (£, V^) on the right of (|3.33p comes from the expansions 
(l2~T2l) . (|375|) and 

— +«7rV = — — + -2ni — - it 2 
dz J h 2 dz 2 h dz 

in other words, L\ (£, Vf ) appears as a coefficient on h^ 1 in the decomposition of A y + (d z + in) 2 in the 
stretched curvilinear coordinates £. Second, 

«4 (£; /S) = «4 (f ; /?) + A± (/3) (2tt |£| 3 ) ft = O (|£p 3 ) (3.34) 

while, according to the rearrangement procedure mentioned above, the main asymptotic term in (|3.9p is 
detached from the right-hand side of (|3 . 33[) because the expression 

A± W) Ci (£, V £ ) = h 3 A ± ((3) d (n, s, z, d n , d s ,8 z ) -^3 (3.35) 

has too slow decay O M£|~ 3 ^ at infinity and, hence, putting (|3.35f> into (|3.33|) would lead to the in- 
sufficient decay rate of the boundary layer term. The transmission of certain unsuitable constituents 
from one limit problem to the other limit problem and the preservation of the behavior of lower order 
asymptotic terms as r — > + and |£| — > +oo implies the absence of the rearrangement procedure. Recall 
that, indeed, the expression (|3.35[) with the cut-off function x is a part of the right-hand side in (|3.11|) . 
and notice that the detachment made in (|3.33|) helps crucially to derive necessary estimates in §3J 

Similarly to Rcmark[8]onc, based on a general result in [18] (see also [19l §3.5, §6.4]), may conclude 
the existence of a unique decaying solution of the problem (|3.33|) , (|3.32[) and the relation 

l 4(f;/3) = 0(|er 3 ), |€| -> oo. (3.36) 

We emphasize that, due to the factor h 2 instead of h in w±, the decay rate (|3.36|) damps down the 
influence of w± on U± (ci. a calculation in (|3.9p ). 

3.6 Simple eigenvalues 

If 77 € [0, 2tt) and r\ ^ n, the eigenvalues (77) = M\ + (77 — it ± 7r) 2 (see (|2.5|l ) are simple and the 
corresponding eigenfunctions U± are still given by (|2.8[) . The asymptotic structures remain the same as 
for (|3.1[) but loose dependence on the parameter (3. In particular, 

w±(0=d n V(O')W(0 (3.37) 

(cf. (|3.8p with a± = 5 a .±) and t/j_ satisfies the equation 

- (A y + (0 Z + 7 ?? )) 2 U' ± (x) - A° ± (77) U' ± (x) = A° ± (77) Ui (x) + 

(3.38) 

+ (A, + (d z + 7T?)) 2 ( x (x) , .xetn , (3.39) 
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supplied with the Dirichlet conditions (|3.12j) and the periodicity conditions. Since the eigenvalue Af]_ (77) 
is simple, only one compatibility condition must be verified, and repeating the calculation (|3.28[) . (|3.29[) 
brings the equalities 

A' ± ( V )=V: =P e \d n V l {0')\ 2 . (3.40) 

One readily sees that formulas (|3.2j) . (|2.5j) . ()3.40j) with 77 = tt + j3h 3 bring about an expansion for the 
eigenvalues A^_ (tt + /3/i 3 ) which differs from the expansion obtained in the previous section. This lack 
of coincidence originates in ignoring the second compatibility condition, namely the norm of the inverse 
operator in w h restricted onto a subspace of co-dimension 1 grows when 77 — > tt and the eigenvalues 
A(j_ (77) approach one the other. 

In $4] the most attention is paid for an appropriate estimate of the remainder k\ (tt + /3/i 3 ) in a 
sufficiently wide range of the deviation parameter (3. 

4 Justification of the asymptotic expansion 

4.1 The operator formulation of the cell problem. 

To estimate the asymptotic remaiders A± (77) in formulas (|2.18|) and (|2.13|) . we employ the following fact 
which is known as "Lemma on almost eigenvalues and eigenvectors" and can be found in, e. g., [31( 13] 
with much more general formulation. 

Lemma 9 Let H be an Hilbert space and K be a compact self-adjoint positive operator in H. If y £ H 
and (f € meet the conditions 

||y;H||=l, \\Ky-<py;H\\=8e(p,<p), (4.1) 

then the segment [ip — 5, ip + 5] C M+ contains an eigenvalue ip of the operator K. 

The space Hp er (zu h ;j h ) equipped with the scalar product 

(u, v) v = ((V, + ir,e 3 ) U, (V, + Z7 7 e 3 ) V)„ h (4.2) 

is denoted by TC(r]). Here 7/ £ [0, 2ir) while the Friedrichs inequality (cf. the middle part of (|2.2ip ) 
provides the positiveness of the Hermitian form (|4. 2|) . 

By a simple argument, the operator K, (77), determined by the identity 

(AC (77) U, V) v = (U, V) mh , U, V £ H (t?) , (4.3) 

is compact, self-adjoint and positive. Owing to [U Thm. 9.2.1], the spectrum of this operator consists 
of the essential spectrum {0} and the discrete spectrum 

$ (V) >^(ri)>-> $ fa) > ... - 0+. (4.4) 
Comparing ()4.2[) . (|4. 3[) with (|2.16|) . we observe the relationship 

A h p (v)=ri(vr 1 (4.5) 
between entries in the eigenvalue sequences (|2 . 1 7|) and (|4. 5[) . 

4.2 Approximation solutions for the spectral problem. 

Let us consider the most interesting case (|3.1|) . We shorten the notation as follows: 

Hp = H (tt ± ph 3 ) , Kp = K (tt ± ph 3 ) , ( , ) n+0h3 =(,)«). 
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Furthermore, we set 

<P± = li 1 ,y± = \\Y ± ;H li \\- 1 Y k , (4.6) 

where 

l± = A + /i 3 A' ± 09) , (4.7) 

y± (x) = X h (x) E/£ (x; (3) + (1 - X h (re)) (n9 n C/£ (0; /3) + U± (a; /?)) + (4.8) 
+ h x (x) (wi (£; 0) + ft«4 (£; 0) + /i 3 X h (ar) (x; 0)) . 



Notice that the dependence on h and (3 is not indicated in (|4.6[> . Addenda on the right of (|4.7|l have 
been determined in (|3.2|> . (|3.31|) . However, the function (j4.8fl has still to be specified. First, the regular 
terms U±, U± and the boundary layer w± were constructed in $2]while the coefficient column a in the 
linear combination (|3.4|) had to be an eigenvector of the matrix (|3.30|) , 



a 



± = a h (v, -2tt/3 ± ^V 2 + 4tt 2 /3 2 ) , (4.9) 



a = 



2 (V 2 + 4tt 2 /3 2 t 2tt(3 v / V 2 + 4tt 2 /3 2 ) , 
(a ± ,a ± ) R2 =l,(o ± ,o =F ) Ra =0. 
Second, the cut-off function A^ is determined in l|2.24|) while, according to (|3.7|) . we set 

U± (x; /?) = 4 (/3) U+ (n, a, z) + o± (/?) U" (n, s, z) . (4.10) 
Third, Xh cuts off the regular terms near the cavern 9 h and, thus, due to the relation (see section [3]3) 

w\ (t;f3) = -t 1 d n Ui(0),ZedQ, 

and the boundary condition (|3.32[) for w± (£; /3), the function Y± vanishes on the surface j h and, therefore, 
falls into Hp er (■cu h ; 7 /l ) . We finally mention that Xh smooths down the correction term U'± which gets 
a singularity at O (see Section [3)3). 

Calculating the norm ||U±; W/3||, we obtain 

(U° ± , U° ± l +0h3 = || (V, + i (tt + ^)) Ul-L 2 || 2 = 
= \\(\7 x +nr)U ± ;L 2 f + R° ± = 
7 y 



V^;L 2 (w)|| 2 +^ 2 ||^;L 2 (w)f+i?^ = 



=Mi +7T 2 + , 

|-R±| <c4(l + 4), (4.11) 
4 =||v :c f/£;i 2 (Mnj' 1 ))! + |tt + /?/i 3 | ||t^;L 2 (ro\ra7 fe ) || + 
+ |/3|/i 3 ||^;L 2 <ch% (l + \0\h* 



and 



\\(l-X h )(U°(;(3)-nd n Ul(O;f3)-XJ ± );H \\ < (4.12) 



< c 



2c x /t 



(I 



(V 4 + (h- 2 + |tt + /3/i 3 | 2 ) r 6 ) r 2 dr J < 



< 



ch* (l + |/3|/i 5 ), 
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h" 



I p 



;Hp\\ <ch p (hi ||V^±;i 2 (6)|| + (l + |j9| h 3 ) \\w±; L 2 (0)\\) < (4.13) 



< ch p+ i (1 



P= 1,2, 



\X h U' ± ;H \\ <ch a U ^ ( r_4+ ( 



°) 



\TrPPf3h a \ 2 ) r- 2 



I c x h 

<ch'i (l + \0\h 3 ) 



3|^.-2\ r 2 dr 



< 



(4.14) 



Let us comment on the above calculations. In (|4.11[) and then in ()4. 12[> we applied the explicit formulas 
([3~4]) , (|49|) , (fL~6j) and also the relations (|3~6j) , (|2~24|) . The inequalities ([443]) hold true due to the 
coordinate dilation a; i — ^ and the inclusion w± G -ff 1 (@) inherited from the expansion (|2.15| and the 
relation (|3.36|1 . Finally, the term h 3 U± with Xh was treated by means of the estimates p.26|) taking the 
properties (|2.24p of the cut-off function into account. 
Imposing the restriction 

|0| < h-^Po, (4.15) 



which damps down the parameter (3 in all bounds in the inequalities f|4.11|l — (|4.13p . We emphasize that 
the weaker restriction |/3| < f3\h~^ is sufhcient here (see the last estimate in (|4.11|) ). however in the 
sequel we need (|4.15[) and just this restriction has been imposed in Theorem|4] We observe (|4.8jl — ()4.14l) 
and conclude that, for a small /3o > in (|4.15|) , the following inequality is valid: 

\Y±\n p \\ > ^(Ah+n 2 ). (4.16) 

Moreover, by (|4.7|) and (|3.31[) , under the same condition (|4. 15|) , the numbers (|4.7p are subject to 

l±>\{M 1 +n 2 ). (4.17) 

4.3 Justifying the asymptotic expansions of eigenvalues. 

For the approximate solution (|4.6p . the quantity 5 in (|4.1[) takes the form 



\jc p y± - <p±y±;H \\ = sup (ic fi y ± - <p±y±,v) 



(4.18) 



Ig 1 \\Y±;Hf,\\- sup (Y ± -l ± !CpY ± ,V) 



; (/3) 



(4.19) 



= Ig 1 \\Y ± ;H P \\ 1 Bup\S±(y)\<csuv\S±(y)\, 

where the supremum is calculated over all functions V G Tip such that || V; || = 1 and 

S± (V) = ({W x +i (r) + (3h 3 ) e 3 ) Y ± , (V x + i (tt + /3h 3 )) V x V)^ h - 
-{K Q + h 3 K' ± {p)){Y ± ,V) mh . 

Notice that the Friedrichs and Hardy inequalities (see (|2.21|) and (|3 . 1 6[) with a = 1) provide the estimate 

\\V X V;L 2 {w h )\\ + Wp^ViL 2 (zu h )\\ < c\\V;Hp\\ = c. (4.20) 

We extend the test function V by null onto w and subtract from (|4.19[) the following scalar products: 

((V, + nre^) U° ± , (W x + «re 3 ) V) w - A {U° ± ,V)^ = 0, (4.21) 
h 2 {X7 e w 1 ± ,X7 e ( x V)) zo = 0, (4.22) 
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h 3 (V f «4, V 4 ( X V)) e - h 3 (£ 1 w 1 ± ,xV) e = 0, (4.23) 
((V, + ine 3 ) U' ± , (V a + i^e 3 ) PW) W - A - 

-A' ± 09) (t/°, + 2z/3 ([/£, + «r) W)) ro - (4 . 24) 

-A' ± 09) ((V, + i7re 3 ) X (2^r 3 ) _1 n, (V, + i7re 3 ) = 0. 

The equality (|4.21| is just the integral identity (|1.17| . The function \V is written in the stretched 
curvilinear coordinates £ (see (|2.9|) ). vanishes on 99 and has a compact support; thus (|4.23| and (14.221) 
follow from the equation (|3.33[) and the harmonicity of the function (|3.34p . respectively. Finally, (|4.24D 
is but a consequence of p. lip . (|3.12p ; note that the test function X^V vanishes near the point O where 
U'± has the strong singularity (|3.25|) . 

In the next two sections we estimate terms which are left in (|4.19p after subtracting left-hand sides 
of (|4~2T]) - (f4~24"P and obtain the common bound cbA \V;Hp\. By virtue of (j4~16)l , (|4~T7)) and flUTS), 
Lemma E] delivers an eigenvalue i\j± of the operator K,p such that 

\il>± -l± \ < c^h*. 

Using (|4.5p and (|4.7p , this formula yields 

|A£ (tt + /3h 3 ) - A - h 3 A' ± (/3)| < Ci ,hik h ± (tt + (3h 3 ) (A + h 3 A' ± (J3)) , (4.25) 



A'l (tt + I3h 3 ) (l - (A + h 3 A' ± 09))) < A + h 3 A' ± ((3) . 



Thus, recalling the condition (|4. 1 5[) and choosing ho > such that the factor A^_ {tt + (3h 3 ) on the left 
of (|4.25p is bigger than — we arrive at the estimate 



| A£ (tt + f3h 3 ) - A - h s A' ± 09) | < c A hl 



which proves Theorem 01 



4.4 Discrepancies of the regular terms. 

Proceeding with U± , we have to take into account the scalar product 

h = h 6 (3 2 (U±,V) m0 

and the cut-off function Xh in (|4. 8[) resulting in 

h = ((V a + i (tt + /% 3 ) e 3 ) (1 - X h ) ([/£ - nd„L/£ - U±) , (V x + i (tt + /?/i 3 ) e 3 ) ^) ro0 . 

Other constituents of (|4.19p . involving C/^, are included to either (|4.2ip . or (|4.24[) . Recall that the test 
function V is extended by zero on the whole cell w° . In view of formulas (|3 . 3[) . (|4.9|) and (|4.20p . a bound 
for |ii| looks as follows: 

2 



ch 6 {3 2 \\V;L 2 (w°)\\ < chi (hip}' 



Here we have used the restriction (|4.15p on the parameter f3, which also applies in further calculations. 
In the sequel we skip mentioning this argument. 
By ([311), (EH and I02QJ), we have 

(r 4 +(h- 2 (l + \(3\h 3 ) 2 )r 6 )r 2 dr\ \\V; Hp\\ < ch* . 
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We now consider the terms due to the transportation of Xh from U± to V, namely 

h = h 3 (U' ± V x X h , (V x + i (tt + f3h 3 ) e 3 ) V)^ + h 3 ((V x + * (tt + f3h 3 ) e 3 ) ^, WJ^ 
We obtain 

r 2c T h 



\h\ < ch 3 



~- 2 r 2 dr || V;Hb\\ + 



+ 



2c~h 



r- 4 + (l + \/3\h 3 ) r- 2 )r 2 drh 



< 



< ch 3 (h^h + h^h 2 ^ 2 ) \\V;Hpf < chi. 



Here we used the estimates (|3.26p and (|4.20[) for U± and V, respectively, while the factors h~ 2 and h 2 are 
caused by the differentiation of the cut-off function and the relation cx < h r < 2cx on supp|Va;X/j| 
(see (|2.24[) and compare with (|2.25p ). The list of other remaining terms reads 

i?h 6 ((t4, (V, + i (tt + (3h 3 ) e 3 ) X h V) m0 - ((V, + nre 3 ) U' ± ,X h V) m0 ) 
-h e k' ± {(3) (U' ± ,X h V)^ . 

7 

Estimates for these terms with the bound ch? become evident after applying the inequality 

\\X h V;H p \\ <c\\V;H \\ 

following from f2~2~4"|) and P~2"0|) . 

4.5 Discrepancies of the boundary layer terms. 

First of all, we replace w± by w± since the main asymptotic term, subtracted in (|3.34[) from the boundary 
layer solution w±, has been included into the equation p. lip and, therefore, the expression (|4.24[) . 
Next, the inequality 



h | ({w± + hw 2 ± ) V xX , (V a + i (tt + I3h 3 ) e 3 ) V) w0 + 

+ ((V a + i (tt + (ih 3 ) e s ) + hw 2 ± ) , FV xX ) . | 



< 



< ch 





r 




r 




I ( 




\h 2 






/supplVjcX \ 


h 




h 


H 



\V;Hp\\ < ch 4 



permits for transporting the cut-off function from w± to V. Note that derivatives of x vanish in a 
neighborhood of the point O and the integral over the support of | V^xl converges. Moreover, the relations 
(|3.34j) and (|3.36|) for w± and w±, respectively, can be differentiated in the case SsupplV^xl that 

was used in the inequality. 
Finally, we write 



h | {i(3h 3 w l ± , (d z + i (tt + (3h 3 ) ) X V) „ 1 + h | ( (d z +i(Tt + ph 3 ))wl, i(3h 3 x V) ro0 1 + 
+h 2 | ( (i (tt + f3h 3 ) w\ , (d z +i{-K + [3h 3 )) X V) ) a0 1 + 
+h 2 \((d z +i(ir + (3h 3 )) w 2 ± , i (tt + /3/i 3 ) x^) CT c | < 
< ch 2 (1 + \/3\ h 2 ) (hi {\\w\-L 2 (6)|| + \\w 2 ± ;L 2 (6)|| 2 ||F;^||) + 

+hi (|||e|V^i;L 2 (6)|| + |||f|V € «4;L a (6)||) Hr" 1 ^ 2 || 



(4.26) 



< 
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h I (Ac + h 3 A' ± (/?)) + ft«4, xF) ro0 1 < 

<ch (hi |||e|wi;i 2 (6)|| Wr^V-.L 2 (m°)\\+hh% ||u4;L 2 (9)|| \\V;L 2 < chi 



Here we have made the transform x i— ► £ which brings the factor ft, 2 on the L 2 — norms of w±, |£| V^w^., 
w 2 -, |£| V^u;^ and the factor /i^ on |£| w±. We emphasize that all norms in L 2 (9) figuring in (|4.26ll 
appear to be finite due to the relations (|3.34|) and (|3.36|) . 

The above considerations demonstrate that the inner products involving boundary layer components 
in (|4.19[) . can be changed with the error O ( hi ) for the sum of the following integrals 



r 1 

h 



mh\ I I (w±d z ( X V) - d z wlxV) (1 + nx) dndsdz, (4.27) 

JO JdujJO ^ ' 

(d n w\d n { X V) + (1 + nx)-' 2 d s w 1 ± d s { X V)+d z w 1 ± d z { X V)) (1 + nx) dndsdz, (4.28) 
(d n w 2 ± d n ( X V) + (1 + nxy 2 d s w 2 ± d s ( X V) + dzwldTixV)) (1 + me) dndsdz, 

(4.29) 

where 1 + nx (s) is the Jacobian, the differential operator V x in the curvilinear coordinates takes the 
form 

— ,(l+nx(s)) — ,t- 



JdioJO 



s 9n ' ds dz / 

and d > is chosen such that suppx belongs to the d— neighborhood of the contour du> allowing for the 
curvilinear coordinate system (n,s,z). 

Replacing 1 + nx (s) by 1 in (|4.29|) brings an error which does not exceed 



ch 2 



\V X {x{x)V(x))\n\h- 1 V^w% (0| da: 



< ch 2 h? |||e| W s w 2 ± (£) ;L 2 (6)|| \\V;Hp\\ < chi 



The resultant integral (with nx = in (|4.29p ) turns into the first addendum in (|4.23|) by the transform 
(71,5,2)1 — >£. 

The same procedure works for (I4.27P with an error less than 

ch f n (\wi (0\ IV, ( X (x) V (x))\ + r \h~ l V ( wi (£)| r' 1 \x (x) V (x)\) dx 

V;H P \\ <chi 



< ch (hi || |£| w\ (0 ; L 2 (9)|| + h^ V&l (£) ; L 2 (9) 
while the norms in L 2 (9) still stay finite in view of 

iei|^(£)ho(ier 2 ), iei 2 |v^(o| = o(iei 

(cf. p.34p ). The resultant integral becomes 

^ f UM-W d i^) = -2*ih* I xV d i±dt. (4.30) 

In (|4.28p we substitute 1 + h£ix (C) and 1 — h^\x(0') for 1 + nx (s) and (1 + nx {s))~ , respectively. 
The concomitant error again gets order hi but, in addition to the expression on the left of (|4.22p . we 
obtain the integral 

h X[0 ) / £1 ~r7Z 7JT 7T C WT- + -JTZ TTT- <% (4.31) 



* , ,s f ( ~i dwi d 2 wl\ 



22 



Since w± is a harmonics, the integrals (|4.3ip and (|4.30j) . according to Q3.33p . form the second term on 
the left of (14T231 . 

We have verified the fact which had been announced in the end of Section[4j3. Our proof of Theorem 
[4] is now completed. 
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